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Abstract 

Let X be a squared Bessel process. Following a Feynman-Kac approach, the Laplace 
transforms of joint laws of (U, Xg ds) are studied where R y is the first hitting time 
of y by X and U is a random variable measurable with respect to the history of X until 
R y . A subset of these results are then used to solve the associated small ball problems 

for Xg ds and determine a Chung's law of iterated logarithm. (^Jq V X§ ds^j is also 
considered as a purely discontinuous increasing Markov process and its infinitesimal 
generator is found. The findings are then used to price a class of exotic derivatives on 
interest rates and determine the asymptotics for the prices of some put options that 
are only slightly in-the-money. 
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1 Introduction 

Let X be a squared Bessel process which is the unique strong solution to 

dX t = 2(u + l)dt + 2^/X~ t dB t , 
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where v > — 1 is a real constant and B is a standard Brownian motion. Letting 6 = 2(z/ + 1), 
X is called a (^-dimensional squared Bessel process. We will denote such a process with 
X = x by BESQ s (x) and 5 and z/ will be related by 5 = 2{y + 1) throughout the text. In 
this paper we are interested in the integral functional 



where p > — 1 and R y := inf{t > : X t = y} for y G [0, oo) (In the sequel, we will write R y 
only if we need to specify the dimension to avoid ambiguity.) and X is BESQ 5 (x). 

Squared Bessel processes have found wide applications especially in Finance Theory, 
see Chapter 6 in [7j for a recent account. They can, e.g., be used to model interest rates 
in a Cox-Ingersoll-Ross framework. In the above setting, if X p models the spot interest 
rates, then exp (E^ ) refers to the cumulative interest until the spot rate hits the barrier 
y p . As such, this random variable is related to certain exotic options on interest rates (see 
[I] for some formulae regarding barrier options in a similar framework). Bessel processes 
also appear often in the study of financial bubbles since l/y/X is the prime example of a 
continuous (strict) local martingale when X is a BESQ(3) (see, e.g., [13], [15] and [16] 
for how strict local martingales, and in particular Bessel processes, appear in mathematical 
studies of bubbles). 

In Section 2 we will determine the joint law of (U, E^ ) by martingale methods, where 
U is a random variable measurable with respect to the evolution of X until R y . In particular 
we will obtain the joint distributions of (R y , E* ) and (max t < Ry X t , E* ). As a by-product 

of our findings, if -Mr = i , we have a remarkable characterisation of the conditional law of 

° ' p+1 2 ' 

^pxy wnen x > V (resp. x < y)) given that the maximum (resp. minimum) of X at R y is 
below (resp. above) a fixed level in terms of the first hitting distributions of a 3-dimensional 
Bessel process. 

We will use the results of Section 2 in order to study small ball probabilities for E^ 
in Section 3. Solving the small ball problem for E^ x amounts to finding the asymptotic 
behaviour of — logProb(Ep < e) as e — > 0. We will then use this asymptotic form to 
determine a law of iterated logarithm for (Ep 0y ) y >o as y — > oo. 

Section 4 will analyse (Ep j/ )i/>o as a Markov process indexed by y and compute its 
infinitesimal generator when v > 0. We will also consider the process Z 6 which is obtained 
via a 'time reversal' from (Ep y ) y >o- More precisely, we will find the generator of Z s defined 



where L x := sup{t > : X t = x}. In particular, we will obtain that Z 4 is identical in law to 
an increasing family of hitting times of a linear Brownian motion. 

Finally, in Section 5 we will apply our findings to the pricing of some exotic derivatives 
on interest rates. The small ball probabilities will be used to find asymptotic behaviour of 
some put options with small strikes, the options that are only slightly in-the-money. 




(1.1) 



by 
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2 Preliminaries 



Let (O, J 7 , (J-f)t>o, P) be a stochastic base where J 7 is completed with the F-null sets. Let X 
be an M + -valued semimartingale which is the unique strong solution to 

dX t = 2(u + l)dt + 2^/X~ t dB t , (2.1) 

where v > — 1 is a real constant, B is a standard Brownian motion. 

X = (flyJ 7 , (J r t), (X t ), (Q 5 x ) x m+) is a Markov process with values in R + . Here Q s x is the 
law of X starting at x. It is well-known (see Chapter IX of [H]) that for v > the set {0} is 
polar, otherwise it is reached a.s.. Moreover, the process is transient for v > and recurrent 
otherwise. We will denote the first hitting time of for X with R. The scale function, s u , 
for BESQ S is given by 

s u (x) = -x~ u for v > 0, s°(x) = log x, s u (x) = x~ v for v 6 [-1,0). 

We refer the reader to [H] and [5] for a comprehensive study of Bessel processes and relevant 
bibliography. 

In subsequent computations we will follow a Feynman-Kac type approach as in, e.g., [6]. 

Lemma 2.1 Let p > — 1,A > and suppose that u G C 2 , solves the following ordinary 
differential equation (ODE): 

x 2 y" + xy' - y[u 2 + \x 2{p+1) ] = 0, (2.2) 
and is strictly positive on (0, oo). Then, (M^ R ) t > is a local martingale where 

Mt ] := u{^X t )XP exp (-± jT^cfa^ . 

In particular, 

l [t<R] dMi U) = H<nM u) " ^) dB, 

Proof. If we let w(x) :— u(y/x)x~% , it is easily seen that w solves 

2xw" + 2{v + l)w' - ^x p w = 0. (2.3) 

Thus, on [t < R] 

dM (u) = 2Ml u) -{X t )^X t dB t = Ml u) ( - 4=] dB t . 

■ 

The solutions to (12. 2p can easily be determined via the modified Bessel functions, I a and 
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K a , of the first and second kind. Recall that (see, e.g., Section 3.7 in \T7\)) la (resp. K a ) is 
an increasing (resp. decreasing) solution to 

x 2 y" + xy' - (x 2 + a 2 )y = 0. (2.4) 

Moreover, they have the following series representation: 

I a {x) = £ , , 1 , ( X -Y m+a K a {x) = ILM^IM. (2 .5) 

KJ ^ m\T(m + a + 1 ) \2J KJ 2 sin (an) K J 

There exist also integral representations (see p. 172 of [17]) as follows: 

I a (x) = [\l-t 2 ) a -h-^dt, CO- 1 -; (2.6) 

r («+ IK (I) J -I 2 

r (i) (i x ) a r™ , i 

K a (x) = Vri 2 ! / e-^-ir**, «>-i (2.7) 



Then, it is easy to check that the solutions to ( 12. 2 p is of the form 

dKj^ ( -^—V\x p+1 ] + C 2 Iwi (—^—V\x p+1 
p+i \p + 1 / P+i \p + 1 

where Ci and C 2 are arbitrary constants. 



{21 



We now return to determining the joint law of (U, ) for arbitrary positive J^Ry- 

s 



measurable random variables U, where is as defined in (11.1)1 . We will analyse the cases 



of negative and positive v separately. 
2.1 The case v < 

Theorem 2.1 Suppose that v < 0, p > — 1 and /e£ Mo(^) := if — - (^^\^\x p+1 ^j . Then, 
(M t "^) t >o a bounded martingale with 

l [t<R] dMt 0) = h<nM U0) ( 4#T - -^1 dB *- 



Proof. In order to show the boundedness property it suffices to show that 

lim u (\/x)x~2 < oo, lim u (\/x)x~2 < oo. 

x—tO x—>oo 

To see that the first limit is finite, apply a change of variable u = xt in (12. 7p to get for a > 

p fl\ 1 poo 

x a K a (x) = ^ 2> 2 ° / e- u (u 2 - x 2 ) a -^ du, (2.9) 



r(a + |) 
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which converges to a finite limit as x — > 0. This implies that lim^o uo(y/x)x 2 < 00. 

In order to find the limit at infinity, we will make another use of the integral representation 
in (EH). Note that, by substituting t = 1 + (j22) turns into 



1 + - ) d«, 



.2a;7 r(a + i)7 

due to which we conclude that 

lim uo(\/x)x~z < 00. 

In view of Lemma [2. II we now obtain that {M^°J) t >o is a bounded martingale. 
It is well-known that (see, e.g. Section 2.8 in [12J) for v > — 1 

x~% K v {y\x) 



(2.10) 



Q 



exp 



2 R y 



y %K V {^/Xyy 
x~%L v (y\x) 



y < x; 
y > x. 



(2.11) 
(2.12) 



y siuiyfiv)' 

Note that the above formulas are still valid when v > 0. 

Since R < 00, a.s. when v < 0, the following is a straightforward corollary to the theorem 
above for < and y < x. 

Corollary 2.1 Letuo be the function defined in Theorem \2. II and suppose that v < 0, p > — 1 
and y < x. If U is J- 'R y -measurable, then for r > 0, 



Qt 



u o(Vx) ( x 



P 5 X U ° [exp {-rU)\ 



where P°' u ° is defined by 



Ur 



dX t 



u Q WX t ) 

for some P x ,u ° -Brownian motion (3. 



M^°\ Moreover, under P x ,u ° , X solves 



t < R, 



+ 1 dt + 2 v /X t d(3 t 



(2.13) 



Remark 1 Note that u (x)x v is decreasing. This follows from the recurrence relation (see 
Section 3. 7 in \Tfl) 



CY 

K'Jx) = --K a (x) - K a ^{x) 



Remark 2 By taking U = Corollary \2.1\ yields the law of S£ for y < x. Comparing 

J$ where 5* = 2(^ + 1), y* = ^ 



this Laplace transform with 112. shows that £ 
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.,i+p 



and Rtl is the first hitting time of y* for some BESQ 5 * 



Moreover, this equality in 



/aw would be valid for y > x, too. These facts also follow from the time-change result given 
in Proposition XL 1.1 in [Lffl . 
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In particular if U = R y , we obtain the following in view of the well-known formulas for the 
Laplace transform of the hitting times of one- dimensional diffusions: 

Corollary 2.2 Let u be the function defined in Theorem \2.1\ and suppose that v < and 
y < x. Then for r > 0, 



Qt 



Uo(y/x) ( X\ 2 $(x 



MVv) \yJ ®(y)' 

where $ is a continuous and decreasing solution of 

2xv " + 2 ( <(V^vg + 1 ) v > = rv _ ( 2.i4) 



Corollary 12. II also allows us to compute the law of ^p X y on the event that a certain boundary 
is yet to be reached. 

Corollary 2.3 Letuo be the function defined in Theorem \2. 1\ and suppose that v < 0, p > — 1 
and y < x. Then, a scale function of the diffusion defined by 112.1 Sty is 



s (x) = / — 2 dy, x > 0. (2.15) 



Thus, for any a > x 



Ql 



^8 



uo(y/x) ( x\ 2 s (x) - s (a) 



PROOF. The representation of the scale function is due to the well-known formulas for the 
solutions of SDEs, see, e.g., Exercise VII. 3. 20 in [T3]. Note that the function is well-defined 
at x = 0. Indeed, it follows from (12. 9p that \im y ^oUo(y/y)y~2 > 0. Thus, lim^o y u o{-\/y) > 
Cy 1+U for some C > 0. Since y^") is integrable for v e [—1,0), the claim holds. 

The second assertion follows from Corollary 12. II after taking r = 1 and U = logl[R Q> # y ] 
via the defining property of scale functions, see Definition VII. 3. 3 in 



Remark 3 The above result in fact gives us the joint law of (max t <# a X t , ). Indeed, 
for any a > x 

[maxX[ < a] = [R a > R y ]. 

Since K\(x) = \p^,e~ x we have more explicit formulas when = — |. 
Corollary 2.4 Suppose that = — ~ and p > —1. Then, fory < x we have the following: 
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i) 



q: 



exp I —rR y — — S 



exp 



x — y 
2v 



§(x) 



where § is a continuous and decreasing solution of 

2xv" + 2 (z/ + 1 + v/aV") u' = ru. 

ii) The function s~o is, up to an affine transformation, given by exp ^— v^A^—j- 
For a > x ; 



(2.16) 



l[fl a >R !/ ] exp 



2 P.a.V 



exp y/X 



x — y 
2u 



sinh (-yf\ a - 
sinh f-x/A^^" 



exp (~VA^) -exp(-7A^ 
exp (- VA^) - exp (-VA^ 

(2.17) 



Note that the expression in (12.171) yields 

Qi 







exp (—Zfp^yJ 


Ra > Ry 



sinh (-V^ 2 -^ 



%inh (-VA^^ 



using the scale function of X under Q 5 X . Comparing this with (I2.12p for v = 1/2 gives the 
following 

Corollary 2.5 Suppose that = — | and p > — 1. Then, for y < x < a we have that 
the law of ^ S pxy conditioned on the event [R a > R y ] is that of the first hitting time of 
[aT v — y~ v ) 2 /4z/ 2 by a 3- dimensional squared Bessel process started at (a~ u — x~ u ) 2 j ' kv 2 . 



Note that, since lim a _ 



l-x a 



— log a;, when y > 0, we obtain that the above conditional 



laws converge as v — > (and, thus, as p — > —1) to that of the first hitting time of (log y/a — 
log y/y) 2 by a 3-dimensional squared Bessel process started at (log\/a — logy^) 2 . This can 
be viewed as the analogous statement of the above corollary when v — and p = — 1. 

Next we look at the case when y > x > 0. Observe that the function u% as defined in the 
theorem below is still well defined and finite at x = in view of, e.g., the series representation 
of I a in 023]). 

Theorem 2.2 Suppose either that p > and —l<u<0or that v — — 1 and p > 0. Lei 

t*i (x) := ( ^y£ p+1 j. Then, (M^^)t>o is a bounded martingale with 



l [t<Ry] dMl 



(«l) 



«1 



'A* 



/or any u > . 
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Proof. Note that [R < R y ] has a positive probability. Thus, we have to pay attention 
to the behaviour of w(x) = Ui(\/x)x~2 at x — 0. Observe that under our assumptions, 
-fpj > —1, thus it follows from the series representation of I a that w(0) > and is finite 
since I a (l) < oo for any a. Next, we will show that w has an absolutely continuous derivative 
over [0, oo). Using the recurrence relations (see Section 3.7 in [TT] ) 

I a -i(x) +I a +i(x) = 2I' a (x), and 

X 

- (I a -i(x) - I a+ i(x)) = al a (x), 

we obtain that 

I' a (x) = I a+1 (x) + -I a (x). (2.18) 
Using this identity it follows from direct calculations that 



VA p-v-i ( y/\ 



p+1 



where 7 = 1 + Since the leading term of I 1 (j^x^ j as x — > is x + 2 + , we see that 

lim^o w'(x) = when p > 0. Therefore, we obtain immediately from the ODE ( 12. 3 j) that 
when p > lim^o xw"(x) = for z/ £ [—1,0). On the other hand, when p = and v > — 1, 



tu'(x) A 1 J ll+ i(vAx) 

M^) ~ 2^ 4(7x5 

by another application of (12.181) . Thus, in view of the series representation of I a 

lim^H = ^hm J ^ (x) 



x->o w(x) 2 x->-o xl v (x) 

Ar(z/+1) A 



4T(z/ + 2) 4(^ + 1) 



Consequently, lim^oo 2(u+i)w' — = since w(x) > for all x > 0. Again, it follows from 
the ODE ( 12 .3p that lim^o xw"(x) = 0. However, this condition implies that f Q w"{y) dy 
exists and is finite. Since this integral equals w'(x) — w'(0) for any x £ [0, 00), we conclude 
that w' is absolutely continuous on [0, 00) and w'(x) = w'(0) + Jq w"(y) dy for any x £ [0, 00). 
Then, in view of Problem 3.7.3 in [8] we immediately deduce that 

w{X t ) = w{X )+ [ 2w'{X s )y/x~ s dB s + [ {2w'(X s )(u + l) + 2w"{X s )X s }ds 
Jo Jo 

= [ 2w'(X s )y/X~ s dB s + ^ [ X%w(X s )ds. 
Jo * Jo 

A simple application of integration by parts formula now shows that M*" 1 ) is a martingale 
with the claimed representation. ■ 



Corollary 2.6 Letu\ be the function defined in Theorem \2. 2\ and suppose that the hypotheses 
therein hold. Then, we have the following for all x < y: 

i) If r > and U is positive and F^-measurable, 



Q 



cxp ( -rU - 2 S lx l2/ 



M 1 (y / x) / X 



where P^' Ul is defined by = Moreover, under P£' Ul , X solves 



dXt = 2 f v^/X^X, + l \ dt + 2 ^Y td(5t 



t < R 



for some P^' U1 -Brownian motion [3. 
ii) For all r > 



Q 



MVv) \yj 

where ^ is a continuous and increasing solution of 

2xv" + 2[ U ' l{v ^^ + l)v' = rv. 



Ml(V^) 

raj t4 sca/e function of the diffusion defined in Ii2.19\) is given by 

r i 

si(x)= / — j-—dy, x>0. 

Thus, for any < a < x 

Ui(y/x) ( x\ 2 s~i(x) — s~i(a) 



Q 



A, 



M i(v^) s i(l/) - s ~i( a ) 



(2.19) 



(2.20) 



(2.21) 



Proof. The proof follows the same lines as in the proofs of analogous results for y < x, 
hence, is omitted. ■ 



Again, since I_i(x) — y ^ cosh(x) we have 
Corollary 2.7 p > and = — ~. Then, for y > x we have the following: 



9 



where ^ is a continuous and increasing solution of 



2xv" + 2 ( v + 1 + \f\x' v tanh 



a/A 
2i/ 



i " I I tf = ru. 



(2.22) 



27ie function §\ is, up to an affine transformation, given by tanh ^— \/A^-J. 
mj For < a < x ; 



<3* 



sinh (-Va^:) 



sinh ( — v^A 2 — : 



2;/ 



(2.23) 



PROOF. Only part iii) needs proof. Note that 



Qt 



A. 



l[i? n >R y ] exp I --S P)a>)B 

On the other hand, 

cosh ( — — ^a;~ 
2v 



cosh ( -^x~ v ) tanhf-^aT" 



tanh ( -^a~ u 

2v 



cosh(-^r-) coshf-^y-- 



tanh ( -^a~ u 

2v 



tanh I -^ x - 



2v 



tanh I — — 



2z/ 



sinh ( cosh f -^a"") - cosh ( -^x~ v ) sinh ( -*£a 



cosh ( -^a" 1 ' 

sinli ( -^Jx-^j cosh + cosh (-^x'A sinh (^a _i/ 



cosh ( -^a~ u 



sinh ( ^2 ( -" - 
cosh (-^a- y 



which yields the claimed representation. 



Note that in fact we do not need to assume p > for part iii) of the above result to hold 
since X is never before R a for < a < x and also that (I2.23j) and (12 . 1 7[) are the same. 
Thus, 
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Corollary 2.8 Suppose that = — | and p > — 1. Then, for y > x > a we have that 
the law of conditioned on the event [R a > R y ] is that of the first hitting time of 

{a~ v — y~ v ) 2 jAv 2 by a 3- dimensional squared Bessel process started at (a~ u — x~ v ) 2 / Av 2 . 

We end this section with a scaling property which will be useful in the subsequent section. It 
is a direct consequence of the scaling property of BESQ S applied to the definition of S^o,,,- 

Proposition 2.1 Suppose that p > 0,u < 0. Then, we have the following identity in law 
for any y > : 

nP +1 y s — y s 



2.2 The case v > 

Recall that when v > the point is polar for X. Thus, one can prove without any difficulty 
that M {ui \ where Ui is as defined in Theorem I2.2[ is a martingale stopped at R y . 



Theorem 2.3 Suppose that p > —1, v > 0. Let u\ be the function defined in Theorem \2.2\ 

Then, (M^^) t >o is a bounded martingale. 

Recall that BESQ 5 is transient when v > 0, thus Q x (R y < oo) = 1 whenever y > x. 
Consequently, we can deduce the following. 



Corollary 2.9 Let ui be the function defined in Theorem \2.2\ and suppose that p > —1, v > 
0. Then, we have the following for all x < y: 

i) lfr>0 and U is positive and J- '^-measurable, 

X 



Qx 



cxp 



" l(V5) ^ r V'[exp(-rf/)] 



where Pz.' ui is defined by -f^- 
ii) For all r > 



M^ l] . Moreover, under P^ , X satisfies MTUX) . 



Q 



exp 



y 2 p> x >y 



U\(y/x) f x\ 2 ty(x) 



MVv) \yj *(y)' 

where $ is a continuous and increasing solution of A2.20\) . 
Hi) For any a < x 



Qi 



— — 

ui{\fx) ( x\ 2 s~i{x) — s~i(a) 



where si is as defined in K2.2V) . 
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Analogous to Proposition [2J] we have the following scaling property. 

Proposition 2.2 Suppose that p > — 1, v > 0. Then, we have the following identity in law 
for any y > 0: 



We now return to the case y < x. 



Proposition 2.3 Let uq be the function defined in Theorem \2.1\ and suppose that v > 0,p > 
— 1. Then, we have the following for all x > y: 

i) If r > an d U is positive and T 'r -measurable, 

A 



Ox 



exp 



where P^ u ° is defined by 

Ry. 

ii) For all r > 



M ( x 0) . Moreover, under P*' u ° , X satisfies (OTEP until 



Qa 



exp ( —rR y 



A. 



Uo(Vx) ( x\ 2 



2 p> x >y 

where $ is a continuous and decreasing solution of \2.1J$ 
Hi) For any a > x 



MVv) \yj ®(y) 



Q 



1 [Ra>Ry] eX P 



2 P' x .2/ 



U 

u (Vx) ( x\ 2 s (x) - s (a) 



uoiy/y) \yJ s o(y) - s Q (a) 



where sq is as defined in $2.15\) . 



Proof. We will only prove i). Note that (AfJJ$) t>o is no longer a uniformly integrable 
martingale since u now explodes at x — 0. Nevertheless, (M^^ y ) t > is still a uniformly 
integrable martingale. Thus, using the optional stopping theorem, we obtain 

u Q {^/x)x^P^ Uo [exp {-rU)\ = Q s x [exp (-rU) M^ o) l 

= Qi 

+Qi [l [Ry=oo] M^] . 

For v = 0, Ry is finite a.s., hence the claim. In case of v > 0, we still obtain the formula 
since u (oo) = 0, and = oo on the set [R y = oo] for X being transient. ■ 

As before, using the explicit form of Ki, one gets 
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Corollary 2.10 Suppose that = \ and p > — 1. Then, fory < x we have the following: 
i) 

gjL.. ( r, \ !f 



exp I -rR y - ^ & p ^ y 



= — exp 

x v 



x v -y v \ $(ar) 



2v 



rv. 



where $ is a continuous and decreasing solution of 

2xv" + 2 (-is + 1 + y/Xx^ v' = 

ii) The function s is, up to an affine transformation, given by exp [VX^j ■ 
Hi) For a > x 

v\ exp ( v^v 



(2.24) 



Q x 



1 [Ra>Ry] eX P 



-Y 5 

'2 p> x >y 



y 

— exp 

x v 



Va 



x - y 
2v 



exp 



exp 



sinh (VAgg 
^ sinh (VA^) ' 



exp ^V^ 1 
(2.25) 



As in the case with ^ = — | we get 

A 







exp 



2 P>z>3/ 



y» y -u_ a -u sinh (y\^f) 
x» x-» - a'" sinh (7a^ 

1 _ (uy sinh (v^^^) 

i-(ir 8 inh (v^^) 

,sinh(v^^) 



a — y 



sinh 



and hence 

Corollary 2.11 Suppose that = | and p > — 1. Then, for y < x < a we have that 

the law of ^p tX>y conditioned on the event [R a > R y ] is that of the first hitting time of 
[a v — y u ) 2 /Av 2 by a 3- dimensional squared Bessel process started at (a v — x u ) 2 /Av 2 . 



Similarly, since Ii(x) = y^ L sinh(x) we have 

Corollary 2.12 p > — 1 and = \. Then, for y > x we have the following: 

13 



where $ is a continuous and increasing solution of 



2xv" + 2 ( -v + 1 + V\x u coth ( — x 

2v 




rv. 



(2.26) 



ii) The function s\ is, up to an affine transformation, given by coth ^v^Af^J ■ 
Hi) For < a < x, 

sinh (VX^f- 



01 



sinh f 



(2.27) 



Remark 4 Comparing parts i) of Corollary \2.J\ and \2.1ft immediately gives us that, for 
x > y, the distributions of ^p tXtV are different when v has different signs. On the other hand, 
Corollaries \2.5\ and \2. 1 1\ imply that they have the same distribution once they are conditioned 
on the event that the maximum of the underlying squared Bessel process is less than a by 
time R y . Same conclusion holds when x < y. 



3 Small ball problem and Chung's law of iterated log- 
arithm 

The small ball problem (also called small deviations) for a stochastic process Z = (Z t ) te j- 
consists in finding the probability 

F[\\Z\\ <e] as e ->■ 0, 

|| ■ || is a given norm, usually LP or L°° . It is connected to many other questions, such as the 
law of the iterated logarithm of Chung type (Chung's LIL for short), strong limit laws in 
statistics, metric entropy properties of linear operators and several approximation quantities 
for stochastic processes. The determination of the above probability is not feasible other 
than in a very few cases and one is inclined to consider the asymptotic behaviour of 

-logP[||Z|| < e) ase^O. 

The solution to the latter problem is also not available in full generality. However, one can 
get this asymptotic behaviour for Gaussian processes (see, e.g., [11] and [10]) or real- valued 
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Levy processes (see pQ). There is a large amount of literature on small ball probabilities in 
the Gaussian setting and one can consult the survey article [TT] . 

As one can expect from the computations made in the previous section, we will be 
interested in the small ball probabilities for the stochastic process (X t )t>o, and the "norm" 



\p,y 



\Z t \ p dt 



where p G (0, oo) and R y = inf{t > : X t = y}. Observe that the above definition is not a 
real norm unless p > 1, however, as the results in this section does not depend on whether 
|| • \\ Pj y is a true norm, this is not a problem. Our results and proofs are close in nature to 
the results of [9]. 

Interestingly, the small ball probabilities for X under the above norm does not depend 
on its index, u, as seen from the next theorem. 

Theorem 3.1 Let X be a BESQ 5 as defined by / TO) and R y = wf{t > : X t = y}. Then, 
one has, for x > and y > 0, 



lim X-hogQl [exp (-X\\X\\%)] 



V2 



A— i>oo 



]ime p logQ s x [\\X\\ p>y <e] 



£^0 



P+1 
1 

_ 2(p + l) 



x 2 



p+i 

y 2 



p+i 

X 2 



p+1 

y 2 



PROOF. Let w(x) = u(y/x)x~% where u = Uq (resp. u = Ui) for y < x (resp. y > x) and 
Mo and u\ are as defined in Theorems 12.11 and 12.21 respectively. Then, it follows from the 
results of the previous section that 



V2\-hogQ 6 a 



Moreover, when u = u , 



exp 



2 II \\p,y 



Vx 



Vx 



.. \ogw(x) 
hm — — 



lim 

A— >oo 



p+i 

X 2 



\ogU (y/x) 

s/X 
hgK 



lim 

p+1 A^oo 

p+i 
X 2 

hm 



pTT Vp+1 



p + 1 



Vx 

\ogKn_ (A) 



p+i 

x 2 



p+i 



p+1 A^oo 



A 



However, using the integral representation (I2.10p in the proof of Theorem 12.11 we obtain 
that for any a > — |, 

\ogK a (x) = C - x - - \ogx + log / e~ u u a ~^ (l + -^-) du. 

2 Jq \ 2x/ 
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log A- m (A) 

Thus, liniA-s>oo f 2 - — = — 1 since 



\v\ 1 / u 

Inn / e"\p+ 1 2 (1 + — ) du 



e u up+ 1 2 du e (0, oo). 



due to x 2 being integrable on [0, a] for any a < oo. This shows that when y < x 

V2 



lim A-hogQf [exp (-A||X||£j] 



A— s>oo 



P+ 1 



p+i p+i 

x 2 — y 2 



log J v (x) 

In order to show the above limit when y > x, it suffices to show that lim^oo — — = 1. 

First suppose v > 0. Using the integral representation (12. 6p we can obtain 



4 (ac 



(x\ a r l 

{2 ( / (i-t 2 )«-i e -^dt + 



r(« + i)r(|) 



r(a + i)r(|) 



7/ i I 

1 ) a - 2 M Q ~ 2 e-"dK, 

2a/ 



where we used the change of variable t = — 1 + - to obtain the second integral above. Thus, 

log I_a_ (x) 
l im 6 v 7 = 1 



x— >oo x 
_ 1 

lim 



I logx + log x y/xl_^_ (x) 



x 



lim 

x— »oo 

= 0, 



log (er x 2~^x^ + ^ f*(l -t 2 )^~h~ xt dt + V2 - ^)^~^u^~^e- u du 



x 



where the last equivalence follows from that lim^oo = 0. Moreover, the last equality 
holds since 

lim e 2 p+ix^ 1 2 / (l-r)p+i 2 e~ dt + V2 / (1 )p+ 1 2 uv+ 1 2 e~ u du 

x->oc J Q J Q 2X 

POO 

= y/2 u^~^e~ u du e (0, oo). 
Jo 

Now, we turn to prove the limit for v < 0. Observe that the since K a (x) converges to 
and I a {x) diverges to oo as x — > oo, the limits we have obtained so far imply for v > 

K'^(x) r_«_ 

lim p+1 . . = lim -£±i- = 1, 

x-±oo K-H—[x) x-hx Ux) 
p+i y ' y ' 

i.e. for any e > there exists an x* such that whenever x > x* we have 

liT, (a;) -AT_^(x)| < eAT_^(x), II' » (x) - I-*-(x)\ < el^(x). (3.1) 

I j+T y ' p+i y n p+l v y ' 1 v y p+i v yl p+i v 7 v 7 
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On the other hand, (12. 5p yields 

logi_>i_ (x) 
l im 6 V ; 



log <! I — * (x) 

p+l 



= lim - 

However, the above limit equals in view of L'Hospital's rule 

lim —EEL— * ^ K 



p+i 



I-^(x) 



2sin(-£ T 7r) 

i" p+ 



»_{X) 
+ 1 v ' 



which equals 1 since for any x > x* 

V 



x 



p+i 



(x) 



i__ii_(a;) — — A_>^ x 

p+l v ' 71" p+1 v 



< £ 



by (13. ip and the fact that K a = K- a for all a > 0. This completes the proof of the first 
assertion of the theorem. 

The second assertion follows by applying de Bruijn's exponential Tauberian theorem (see 



P+l p+l 

W \X 2 -y 2 



Theorem 4.12.9 in |3J) to a = —1 and f3 — 2 ( p +i)2 

Observe that E^ is an increasing process when indexed by y. We will next use the 
above theorem to obtain Chung's LIL for (Ef 



Theorem 3.2 Lei <f)(y) := 



p+i 



log log j/ 



p.O.yy^O- 

. Then, for any v > — 1 and p > one /ias 



E 5 

liminf 3 



0( y ) 2(p + l) 2 ' 
Proof. It follows form Theorem 13. II that 

lim e log Q s [Ej, 0|1 <e] = 



1 



£-5-0 



2(p+l) 2 ' 



thus, for sufficiently small e, 



Ql [Ej >0il < e] < exp f-^ 



where A is a fixed, but arbitrary, constant in (0, 2 ( P +i) i ) • Fix C > 1 an d se ^ Z/n = C n - Next 
choose > so that kC p+1 < K, In view of Propositions 12.11 and 12. 2\ we get for all large n 

r kc p+1 

< exp 



^^loglog^. 



A" 



T log((n + l)logC) 1 = (logC) fccp+i(n + l) 
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which is summable in n. Therefore, by the first Borel-Cantelli lemma, we have that, a.s. for 
large n, ^fffy > k. On the other hand, for y e [y n , y n+ i], 



which shows that 



and thus 



liminf p ' ' v > k, Qg-a.s., 



y^oo <j)(y) 



liminf 



> 



by the arbitrariness of C, K and k. 

We now turn to prove the reverse inequality. First, let's observe that 



liminf p '°f < oo, Qi-a.s.. 

y— >oo y p 



(3.2) 



The above claim follows from a direct application of Fatou's lemma since = y p+1 Ep 01 . 
Next, fix an e > 0, let y n = n n and consider the events 



E n :- 



Xfds < (l + 2e) 



2(p+ l) 2 



0(j/n) 



It follows from the strong Markov property of X that -E n s are independent, and we will now 
see that E n s occur infinitely often due to the second Borel-Cantelli lemma. Indeed, by the 



definition of T,i n „ , we obtain 



Qi{E n ) > Qt 

= Q s o 



> exp 



Ej A , n < (1 + 2e) 

E s pA1 < (1 + 2e) 
l + e 



2(p+l) 2 
1 



<f>{Vr, 



l + 2e 



2(p + l) 2 loglog?/ n _ 
log log y n J > 

J log Vn 



where the second inequality is due to the fact that, for a given e > 0, Q 5 [S*o i < ??] > 

exp ^— (1 + e) 2 (p+i)2 ^ f° r sufficiently small 77 in view of the convergence result of Theorem 

13.11 Since — r — is not summable, it follows from the Borel-Cantelli lemma that £L occurs 
infinitely often. As e was arbitrary this allows us to conclude, a.s., 

S::, x " ds 1 

lim inf —, — ; < 



<f>(Vn) 



2(p + l) 
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Thus, Qq-&.s., 

li m i n f P-°-y" < ii m i n f p- '^- 1 + - i; m ; n f p<°>v»-i ^n-i f /„ ^ 

On the other hand, 

ygl < loglogn^ +1 > 
which converges to as n — > oo. Therefore, in view of (13. 2 p and (I3.3p . we obtain 



hminf-^< - gg-a.s.. 

n^oo 0(y n ) 2Q»+1) 2 



4 Feller property and 'time reversal' 

In the previous section we have proved a law of iterated logarithm for £^ by considering 
it as a process indexed by y. In this section we will see, for v > 0, that it is in fact an 
inhomogeneous Feller process and find its infinitesimal generator. 

First of all, it immediately follows from the strong Markov property of X that (£* y , J z R y ) y >o 
is Markov. Suppose P X}V is the associated semigroup, i.e. P X) yf(a) = Qo[f(^po,y)\^i,o,x = a]- 
Since the increments of (Ep t o, y )y>o are independent, we have for any bounded measurable / 



Px, y f(a) 



f(a + b)Q s x (E d p ^ y 



e db). 



(4.1) 



Let Co denote the class of continuous functions on M + that vanish at and oo. ( 14. ip readily 
implies that when f £ Co, P X yf £ Co as well. Moreover, it follows from Corollary 12.91 



and the observation that u\ is finite at 0, that for each x > the measure Qf.^p^ G d&) 
converges weakly to the Dirac point mass at as y 4- £ since its Laplace transform converges 
to 1. Therefore, lim,,.^ P x . v f(a) = /(a) and consequently (£* „, J r R„,)„>n is Feller 



The form of the infinitesimal generator of 
theorem. 



^p,0,J/> ^Ry 



p,0,3/' J Hyly>0 

y >o will follow from the following 



Theorem 4.1 Suppose that v > 0. Then, for every x > t/iere exists a decreasing function 
ti(x, •) satisfying 



™ e-? b W(x,b)db = ^^l, Ag(0,oo) 
A w(x) 



where w(x) 



x |-^t4x P 2 ). Moreover, 







lim 

24a: 



exp (-AEj Al/ ) 



P|0,!B 



— e 



- Aa 



-A(a+6) 



-Aa\ 



}n(x,db), (4.2) 
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where ir(x, db) := — tt(x, db) for b > 0. In particular, 

/ bir(x, db) < oo. 
Jo 

Proof. In view of the strong Markov property of X, 



exp 



2 pfi'V 



J p,0,x 



exp 



2 P>aM/ 



Thus, it follows from Corollary 12.91 that 



Q 



lim 

ylx 



exp ( |Ep )2/ ) 



e 2' 



y-x 



c -^ lim ^[ex P (-|S^J]-l 

?4z y — X 



e 2 a li m 



ui(x) 



?4x ?/ — X 



l = _ c -#a^.(4.3) 

w(x) 



On the other hand, using integration by parts, we obtain 

A 



Q 



exp 



2 V,x,y 



poo 

1 = / e-^gjp^ G eft) - 1 

/■oo 

/ e^ b Q 6 x (E 6 P;X;y >b)db. 
Jo 



A '' x 
2 



(4.4) 



It is well-known (see Corollary 3.8 in Chap. VII of [14] ) that 



O 5 \T S 

t» x L p,x,y 



z p G(x, z)m(dz), 



where m and G are the associated speed measure and Green's function, respectively. In our 
case, these are given by 

m{dz) = —dz; G(x, z) = -y~ v + (xV z)~ v '. 



Consequently, the above formula yields 



Qx Pp,£,J/] 



,p+i 



x 



P+i 



2( P + l)(p + jy+ 1)' 



which in particular implies that H(x,y,db) := 2 ^p+i_^ p +"i 1 ^ Qt(^p,x,y > fydb is a probability 
measure on [0, oo) for each (x,y). However, ( 14 .3p and ( 14. 4 p imply that H(x,y, •) converges 
weakly as y tends to x to some probability measure, II (x, •) on [0, oo) which satisfied 



e-^ b U(x,db) 



4(p+ u + 1) _ p w'(x) 
A w(x) 



(4.5) 



1 Using the integral representation of I a for a > — I, it is tedious but straightforward to check that this 
representation holds for x — as well by taking the limit as x — > and showing that L(X) < oo for A > 0. 
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Moreover, since Q x {T, s pxy > b) is decreasing in b for each (x,y), the limiting measure II is 
necessarily of the form 

ce (db) + 2(p + u + 1)x~ p tt(x, b) db, 

where Eq is the Dirac point mass at 0, c a nonnegative constant, and W(x, •) is a decreasing 
function for each x. In particular, W(x, oo) = 0. In order to find the constant c, it suffices to 
check the value of the function L at oo. However, using the explicit form of w, 



c = lim 

A— >oo 



4(p + i/ + l) _ p wf(x) 



A 



-x 



lim 2 



p + v + 



1 \p+i 



W[X) 



p+i \p+i j 



since for any a > — \ 



ir a (x) i i - t 2 y-he- xt dt 



Ob 1 ( X ( X 



a— 

X z X 



and that 



1^(1 -t 2 ) a -ke' xt dt 



f\(i-t 2 ) a ~h- xt dt 



< 1. 



Thus, we have shown that 



w (X) 



w(x) yix y — X 



e 2 b 7r(x, b) db. 



Since n is decreasing with ir(x, oo) = 0, we obtain by integrating by parts 



lim 

yix 



exp (-fSj A J 



^p,0,x — a 



e i 



A 



y-x jo 

where n(x,db) = —W(x,db). Finally, note that one necessarily has 

/ bn(x, db) < oo, 
Jo 

since otherwise L, as defined in (14. 5ft . would have been infinite. 



e" 6 -lf 7r(x,d&) 



The above theorem yields that the sequence of measures 



J/ — X 



c?6 ) converges vaguely 



to a finite measure on (0, oo) as y — > x. Thus, for any / G C^-(R + ,1R), i.e. the space of 
continuously differentiable functions with a compact support, we have 



lim Qg [/W,oJ|s;,o, a = q]-/(c 

l4x y — X 



{f(a + b)-f(a)}n(x,db). (4.6) 
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In other words, letting B(R+) denote the bounded Borel functions defined on M + , if we 
define the operator A x : B(R + ) i-> B(R+) by setting 

POO 

A x f(a):= / {f(a + b)-f(a)}n(x,db) for / G C^(R + , R), 
Jo 

then we see that the process (M^) y > defined by 

M f y := f& s pAy ) - [ V A x f(Z s pAx )dx 
Jo 

is a martingale with respect to the filtration (J 7 R y ) y >o whenever / belongs to the domain 
of A x for all x > 0. The form of the infinitesimal generator also reveals the fact that the 
increasing process (££ ^) y >o is purely discontinuous, i.e. there is no interval (a, b) in which 
it is continuous. 

Remark 5 It follows from the fact that (R y ) y >o is left continuous that o,y)y>o is a left- 
continuous process. However, in view of the above Feller property one can obtain a cadlag 
version of it when we augment the filtration with the null sets. Existence of a right- continuous 
version can also be independently verified by observing that 

limQ* [exp (-A£j >0 J] = Ql [exp (-AE^J] ■ 

We will end this section by analysing a specific 'time reversal' example. To this end let 
L x := sup{t > : X t = x} and suppose that v > so that Qq(L x < oo) = 1 for all x > 0. 
We will consider the process Z s defined by 

Z 5 X := I 1 X p s ds Vse[0,l). (4.7) 

J Li- X 

In view of the well-known time reversal results for diffusions, see, e.g., Exercise 1.23 in 
Chap. XI of [H], the law of the process (X Ll _ t ,t < L\) under Qq is identical to that of 
(X t ,t < Rq) under Q\~ 2v '. Recall that Q\~ 2u (Rq < oo) = 1. Thus, we can write 

Z s x = I 1 X X?ds, X = BESQ 2 - 2u {l). (4.8) 
Jo 

Note that the above equality is to be understood in the sense of equality between the laws 
of the processes. Due to the strong Markov property of X we again have that (Z x ) x ^q^ is a 

Markov process with respect to the filtration (J-'aOxep,]-) where T x := o~{X s \ L\^ x < s < L\). 
Observe, more easily in view of (14. 8p . that Z 5 , too, has independent increments rendering 
its Feller property in view of the arguments that led to the Feller property of (£* x ) x >o 
at the beginning of this section. The next theorem will yield the form of the infinitesimal 
generator. Its proof will follow similar lines of the proof of Theorem 14.11 so we will only give 
the details when it differs. 
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Theorem 4.2 Let Z 5 be as defined in ( [^.7p and suppose v G (0, 1] . Then, for every x G [0, 1) 
there exists a decreasing function n(x, •) satisfying 



e z b Tx{x,b)db 



2w'(l - x) 



A w(l — x) ' 

where w(x) := x^K_»_ ( -^-x 2 ^ - ) . Moreover, 
v ' P+l Vp +1 / 



A G (0,oo) 



^ Qg [exp(-AZ;)|Zi = a]-e 



-Aa 



where p[x, > 



—nix, 



y-x 

for b > . In particular, 



{« 



-A(a+6) 



-e~ Aa }p(x,d&), 



(4.9) 



/ bp(x, db) < oo. 
Jo 



Proof. First observe that in view of (14. 8 j) and the aforementioned independent increments 
property 



exp ( ~-Z 5 y 



A 
2' 



Z£ = a 



cxp 



2 P>1~ x,l— y 



where (3 = 2—2v. Therefore, the same arguments in the beginning of the proof of Theorem l4.ll 



yields that the measures — 



(XT , , >b)db 

-x V — x,l — y ' 



H(x, •) on [0, oo) which satisfies 



y-x 



converges vaguely to some probability measure, 
1 w'(l-x) 



e-i b Tl(x,db) = - 



A w(l — x) 



(4.10) 



Moreover, the limiting measure is necessarily of the form 

ce (db) + 7r(x, b) db, 

where Eq is the Dirac point mass at 0, c a nonnegative constant, and tt(x, •) is a de- 
creasing function for each x. Observe that since Qi- X (jt, t \- X i- y ) = oo, the measures 
Qi- x (^ P ,i- x ,i-y> b ) db are no ^. g n ^ e anc j ne ither is their vague limit. On the other hand, we 

y—x ° ' 

can still conclude that c = since, in view of (I2.10p . for any a > — ~ one has 



1 K(x) = _1J_ _ 1 _ a - \ f~ e-«u°+1i (l + £)' 
2x 2 



3 

2 du 



3C -^^"q; 



which converges to as x — > oo. Next, integrating (14. lOf) by parts yields 

w'(l — x) 



1 5 

2 



7r(x, OO) 



(1 - e-t 6 ^ 



tt(x, db) 



w(l — x) 



(4.11) 
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Thus, due to the dominated convergence theorem, taking the limit of the above as A tends 
to yields 

-/ \ ,. w'(l-x) 

7TIX, oo) = — iim ; -. 

A^O W(l - X) 

However, using the representation in (12. 9 j) one has that 



K a (x) ' V 2 J e~ u {u 2 - x 2 ) a ^ du 

for any a > — ~. However, x 2 J°° e~ u (u 2 — x 2 ) a ~^ du converges to as x tends to due 
to monotone convergence theorem if a > 3/2 or dominated convergence theorem otherwise. 
This in turn yields that 



p+i 



p+i„_**(£*"*M 



v W'(X v p+i,. V P+l J VA £±1 
hm — -— = 1 — hm 7 i -x 2 



A^O w(x) 2x 2x A->0 j£ v ( v^EiA p + 1 

and thus 7r(x, oo) = 0. As in the proof of Theorem 14.11 

c ^ a w'(l-x) = ^ Qj [exp(-|z;)|Zf = a]-e-f^ 
w(l — x) yix y — x 

Thus, combining above with ( 14. lip and plugging in the value of n(x, oo) yield 

lim P V 2 W 1 * ^ = / \ e -^ a+b) -e~ a \p (x,db . 

s/ix y - a; Jo 1 J 



Example 4.1 As an application of the above theorem consider the case when v e (0, 1] and 
-^-j- = | . T/ien, i/ie associated n is defined by 

e'^(x,b)db= -^{l-xy- 1 

in view of the explicit form for Ki . Thus, by inverting the above transform, we have 

1 



7t(x, b) = (1 — x) 



u-l 



y/2~nb 

This reveals that the infinitesimal generator, A x of Z 5 is defined by 

f°° 1 
Ax = {l ~ Xr Jo {fia + b) ~ m} 272W db 
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for any f in C 1 with a compact support. Consequently, 

1 

: dbdy 



f*CC (* OO 1 

im-Is-^L {!{z > +i) - ![z * )] w^ 

4- 

is an J 7 -martingale for such f . 



Example 4.2 Observe that although Z s , or (£* ) ; is an increasing process with indepen- 
dent increments, it is not a subordinator (see JM/ for a definition and further properties) 
since the increments are not stationary. However, in the above example, if one takes v — 1, 
then one sees that Z s becomes time homogeneous, i.e. Z s is a subordinator. Observe that 
v = 1 implies p = 1 in this framework. More precisely, 



%x = Z x = / X s ds, 

'Li-* 



where X is BESQ 4 (1), is a subordinator. Moreover, Corollary \2.4\ and yield 



Q 5 o (exp (-A£*)) = exp ^-y^^j 



Thus, = Ts., where T x is the first hitting time of x for a Brownian motion starting at 0. 



5 Applications to finance 



Our aim in this section is to give some examples arising from some financial models and 
discuss how the results from previous sections can be used to obtain prices of certain financial 
products. 

As explained in Introduction the process X is commonly used in the finance literature 
to model interest rates. Suppose the spot interest rate is given by X p where p > — 1 and X 
is BESQ s (x) and consider the following exotic derivative on interest rates which pays one 
unit of a currency at time R y if the accumulated interest is less than k, i.e. < log k. 

This is an example of a digital option and its price, as usual in the Finance Theory, is given 
as an expectation of its discounted payoff: 

D(k- 5, p, x, y) = Q 5 X lp£„ iV <fc] exp (-S* X)W ) 

On the other hand, if one computes the Laplace transform of D(k; S,p,x,y), one obtains 



e ^ k D{k] 8,p,x,y) dk = Q 



POO 

/ e-^exp(-S^) dk 

J yS 

p,x,y 



Qi [exp(-( M + l)S^J] 
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which is at our disposal due to the results from Section [2j In particular, the above identity 
implies 

' e v D(u',8,p,x,y)du = Qi[J: s p ^ y <z], Vz > 0. (5.1) 



Jo 



Moreover, once the function D is determined by inverting the corresponding Laplace trans- 
form, one can also obtain the prices for put options written on the accumulated interest. 
Indeed, for any K > 1 

P(5;x,y,K) := Q s x [(K - exp (E^ J) + exp 

= Q x J l[s| ia:iV <iogfc] dkexp (—E PtX>y ) 

/•K /"log A" 

= / -D(log k] 5,p, x, y) dk = / e u D(u; 8,p, x,y) du. (5.2) 
J\ Jo 



Note that for K = 1, i.e when the option is at-the-money since the cumulative interest 
at i = is defined to be 1, the put option is worthless. We can in fact get how fast the 
option becomes worthless as K approaches to 1. Indeed, comparing ( 15. 21) to (15.11) yields the 
following asymptotics for the option value in view of Theorem 13.11 



limloginogP^x^ff) = -— — — - 

ku 2{p + iy 



p+i 

x 2 — y 



p+i 



(5.3) 



The above expression tells us how small the option price becomes when the option is slightly 
in-the-money, i.e. when K is very close to 1. 

Next, assume y < x and consider another type of a put option on the maximum of the 
short rate with maturity R y and payoff (K — max ( </; s X t ) + for K > x. The price of this 
option equals 



Q 



f K i r r K 

exp (-£j iSe J / l[max t < % x t <a]da = Q 5 X exp(-E^ J / l [Ra>Ry] da 

Jo J . Jx 



u o(Vx) ( x 



So{x) — s{a) 
sb(y) - s Q (a] 



da, 



in view of Corollary 12.31 and Proposition 12.31 where the pair (uq, sq) is computed by setting 
A = 2. 

Finally, if one is interested in pricing an Asian option on the short rate until time R y , it 
suffices to use the identity 



Q 



yS 



exp {-E 5 p x y ) 



Q 



exp (-Sp, x ,J / 
Jo 



and invert the joint Laplace transform of R y and E p obtained in Section [21 
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